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Abstract. Angular fluctuations of stochastic gravita- 
tional wave backgrounds (GWB) produced by extragalac- 
tic astrophysical sources are calculated. The angular prop- 
erties of such backgrounds are determined by the large 
scale structure of Universe (galaxy clustering). The evo- 
lution of star formation rate with redshift is taken into 
account. Fluctuations of the metric strain amplitude as- 
sociated with such noises at angular scales of about one 
degree are found to be of order 5-20% slowly growing to- 
ward smaller angular scales. This feature can be poten- 
tially used to separate astrophysical GWB from cosmo- 
logical ones in future experiments. 
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1. Introduction 

With the advent of detectors of gravitational waves (GW), 
which are currently under construction, a totally new pos- 
sibility to study different astronomical objects opens up 
(Thorne 1988). In addition to "classic" sources of GW like 
merging compact binaries or rapidly rotating hot neutron 
stars etc., which will be studied at frequencies ~ 10—1000 
Hz, there should exist a specific cosmological background 
(noise) covering a very wide frequency band from ~ 10~ 18 
to 10 5 Hz. The cosmological background should bear im- 
prints of the physical processes in the very early Universe 
(Grishchuk 1988, 1997 and references therein). The pri- 
mordial GW background, which originates from vacuum 
fluctuations parametrically amplified by the very expan- 
sion of the Universe, has a power-law spectrum spanning 
a very wide frequency range (Starobinsky 1979, Rubakov 
et al. 1982, Grishchuk 1988, 1997). The prospects for its 
detection appear to be the most favorable at the low fre- 
quency band (10~ 3 — 10 _1 Hz) which will be covered by 
LISA space interferometer (Larson et al. 1999). 

At some frequencies, however, unresolved binary stars 
within our own Galaxy or beyond provide an important 
contribution in the LISA frequency band (Bender and Hils 
1997, Postnov and Prokhorov 1998, Kosenko and Post- 



nov 1998). As shown in these papers, the stochastic sig- 
nal from unresolved merging white dwarf binaries domi- 
nates the LISA sensitivity curve up to ~ 3 x 10 -3 Hz. At 
higher frequencies extragalactic merging white dwarf bi- 
naries contribute at a level roughly 10 times smaller, which 
is below the planned LISA sensitivity at these frequencies. 
This fact allows search for detection of the primordial GW 
backgrounds by LISA (Grishchuk 1997). 

However, some inherent uncertainties (e.g. in the 
galactic rate of binary white dwarf mergings) are present 
in the calculations so the astrophysical backgrounds can 
turn out to be higher than expected. Since astrophysical 
GW backgrounds are considered as an additional noise 
contributing to the intrinsic noise of the detector, as much 
as possible of their properties at all frequencies should be 
known in advance. 

What are specific features of astrophysical GW back- 
grounds? Clearly, those related to the galactic sources 
should follow the distribution of stars inside the Milky 
Way (Hils et al. 1989, Lipunov et al. 1995). As all the de- 
tectors (ground-based or space-born) should rotate with 
respect to the galactic plane, the signal modulation has 
been used as an advantage to detect them (Giazotto et 
al. 1997, Giampieri and Polnarev 1997). As for the back- 
grounds of extragalactic origin, only amplitudes at differ- 
ent frequencies from various sources have been computed 
so far (Kosenko and Postnov 1998, Ferrari et al. 1999a, b). 

The purpose of the present paper is to study angular 
properties of the GW noise produced by extragalactic as- 
trophysical sources at the degree scales. As most of these 
source must reside in galaxies (only a tiny fraction of bina- 
ries or GW-emitting neutron stars is expected to be in the 
intergalactic space) , the GW background should have dis- 
tinctive angular correlation properties reflecting the large 
scale structure (LSS) of the Universe. This is exactly what 
we observe in electromagnetic radiation as fluctuations of, 
for example, IR background observed by COBE (Kashlin- 
sky et al. 1999). 
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2. Imprint of the LSS in the astrophysical GWB 

The angular fluctuations of intensity in any astrophys- 
ical GWB should correlate with those of the projected 
number of galaxies, so first we wish to remind the reader 
the well-developed technique which is used for studies of 
angular properties of galactic counts on the sky. There 
is a lot of specially dedicated literature on this subject 
(see e.g. Peebles 1980, 1993). Before to proceed further, 
let us make a simple estimate. Suppose we have sources 
(galaxies) randomly distributed in space. What is the rms 
fluctuation of the number of galaxies within a small solid 
angle Sil = tt9 2 <C 47r? With the average space density 
of galaxies uq ~ 0.01 Mpc -3 and neglecting cosmological 
effects for a while we would get 
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(here H is the Hubble constant) while galactic counts 
demonstrate a much higher value 5N/N ~ 0.5 for linear 
scales I — 30 Mpc (i.e. 9 ~ 0.01 for the most distant 
galaxies) (Peebles 1993). So clearly we should take into 
account the correlation properties of galaxies and galactic 
clusters. 

To account for the non-Poissonian properties of galac- 
tic distribution in space, the correlation functions are 
used. For Gaussian fluctuations only two-point correla- 
tion function £(r) or its Fourier transform (power spec- 
trum) Pa(fe) would be sufficient. Taking Pa(fc) as derived 
from LSS studies, we then can calculate the 2-dimcnsional 
correlation function C(9) of the projected distribution of 
galaxies on the sky or, equivalently, the 2-dimensional 
power spectrum P2{q)- This is the last quantity that we 
actually need since the rms fluctuation of the energy flux 
per unit logarithmic frequency interval in a given direc- 
tion is directly related to ^2(9) (Kashlinsky et al. 1999 
and below). 

2.1. Correlation functions and power spectrum 

The brightness of a GW background can be characterized 
by the energy flux coming from a given direction within a 
solid angle Afi (Thorne 1988) 
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where x is the two dimensional coordinate across the 
sky. The integration over all sky yields the familiar value 
flcwPcrC 2 , the energy density per unitary logarithmic 
frequency interval in units of the critical energy density 
to close the Universe, which characterizes the isotropic 
stochastic GWB. To study angular properties of the noise 
we shall consider the flux from given direction F(x) = flf 



itself. We shall assume some ideal GW detector with a 
beam-like sensitivity diagram, which is of course far from 
realistic ground-based LIGO-like or spaceborn LISA-like 
interferometers. However, in this paper we will not discuss 
the observability of GW backgrounds. 

The fluctuation in the GW flux arrived at the de- 
tector from a given direction on the sky is <5F(x) = 
F(x) — (F), where (...) denotes ensemble averaging. 
The Fourier transform of the fluctuation is SF(9) = 
1/(2tt) 2 / 8F q cxp(-iq0)d 2 q. 

The projected 2-dimcnsional correlation function rep- 
resents the first non-trivial moment of the probability dis- 
tribution function of <5F(x): C{9) = (5F(x + 9)F(x)). 
The two dimensional power spectrum is by definition 
^2(3) = (\SF q \ 2 ). We shall assume the phases to be ran- 
dom and the distribution of the flux field to be Gaussian 
so that the power spectrum is just the Fourier transform 
of the correlation function: 
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Here Jo is the zero-order cylindrical Bessel function. 

The mean square fluctuation of the flux on the detector 
within a finite solid angle A£l subtended by the angle 1? 
across the sky is zero-lag correlation signal 
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where W is the window function of the detector. For exam- 
ple, for a top-hat beam ((SF) 2 )$ ~ {l/2ii)q 2 P 2 {q)\ q ^ 7T /2V 
and the values of 1/q correspond to fluctuations of angular 
size ~ n/q. 

The GWB flux and its angular properties measured 
in projection on the celestial sphere should reflect 3- 
dimcnsional structure of the Universe and the change of 
GW emission rate with redshift z. The LSS can be taken 
into account by the 3-dimensional correlation function 
£(r) or its 3-dimensional power spectrum P$(k), which for 
isotropic case relates to £(r) through the equation 
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(here j is the zero-order spherical Bessel function). 

The projected correlation function of cosmic GWB 
C{9) is expressed through two-point correlation function 
of the galaxy distribution £( r i2) and the rate of GW emis- 
sion dF/dz via the Limber equation (Limber, 1953): 
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Substituting this equation in the limit of small angles <C 
1 for Friedman- Robertson- Walker metrics into Eq. (Q) one 
gets (Kashlinsky et al. 1999) 
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where dA — d m /(l + z) is the angular distance, d m is 
the metric distance. For degree angular scales of interest 
here the linear approximation of galactic clustering can 
be used P${k, z) ~ Ps(k,0) x *&(z), where $>(z) describes 
the evolution of the clustering. On linear scales ^(z) = 
(1 + z)- 1 if CI = 1 (Peebles, 1980). 

In our analysis we use the 3D spectrum Ps(k) as de- 
rived by Einasto et al. (1999) from a thorough analysis 
of different LSS studies. The mean spectrum of galaxies 
shows a power-law behavior at small and large k with a 
maximum ~ 10 4 [/i~ 3 Mpc 3 ] at k w 5 x 10~ 2 h Mpc -1 . 

The analysis of this formula (Kashlinsky et al. (1999) 
and references therein) shows that the relative fluctuations 
of the flux are SF rms ~ yj H /ck 2 P 3 (k) ~ 5% - 10% and 
weakly dependent on the cosmological model at angular 
scales of order one degree. However, strong dependence 
on redshift of the flux rate requires more accurate calcu- 
lations. 



2.2. GWB flux 

The GW energy emitted by the population of some sources 
in a galaxy at frequency / per unit logarithmic frequency 
interval d In / in the rest-frame of the galaxy can be cal- 
culated through the rate of GW-producing events 1Z (e.g. 
the rate of binary WD mergings or supernova explosions) 
in the galaxy. Under the stationary conditions we have 
(Kosenko and Postnov 1998) 
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where E(f) is the energy which is being carried away by 
gravitational waves from the typical source (e.g. for two 
point masses orbiting each other it is just the orbital bind- 
ing energy). If E(f) oc f a and GW emission is the only 
dissipative mechanism 
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The comoving luminosity at proper frequency /' per 
unit logarithmic frequency interval produced by sources 
in galaxies at redshift z from the redshift interval dz from 
unit solid angle 
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where n(z) is the space density of galaxies, SV(z) is the 
proper volume clement. We assume no new galaxies to 



create since their formation so that n(z) = n G {\ + z) 3 
and n G = 0.013(f2b/0.005)/iio O Mpc -3 is the present-day 
density of galaxies normalized to the amount of baryons 
comprised in stars (in terms of the critical energy den- 



sity p c 
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if /100(km/s/Mpc) is the Hubble con- 



stant). The strong star formation rate evolution is taken 

into account through the evolution of the event rate with 

redshift IZ(z) (Kosenko and Postnov 1998). We use the 

parametrization of Rowan-Robinson (1999) for the star 

formation rate history SFR(z) as derived from optical, 

UV, and IR-observations, and normalized to unity at 

z = 0. The rate of particular events is thus TZ(z) = 
o 

1Z G J SFR(z')G(z — z')dz' where 1Z G is the present-day 

event rate per galaxy, G(z — z') is the redshift dependence 
of the rate after a 5- function- like star formation burst. 
This is important especially for binary white dwarf co- 
alescences because these are delayed typically by ~ 10 9 
years since their formation. For events which relate to 
massive star evolution, like supernova explosions, we have 
G(z — z') w S(z — z') and the change in the rate of events 
with redshift simply follows the star formation rate evo- 
lution n{z) =1l G x SFR(z). 

The proper volume element is (Peebles 1993) 
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with the metric distance d m and dt/dz being the func- 
tions of the cosmological model (Carroll et al. 1992). In 
our calculations we use the standard flat universe without 
cosmological constant and a A-term dominated cosmolog- 
ical model with J7a = 0.7, Sl m = 0.3. 

Combining equations (pp|), (|TT|) , and (|l^) together we 
arrive at 
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so the contribution to the total GW flux from the redshift 
interval dz from unit solid angle which is observed today 
by a detector with band-width df centered at frequency 
/ = /'/(l + z) is 
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where we used the luminosity distance definition d^ = 
d m (l + z). Note that no additional factor (1 + z) account- 
ing for the change of frequency interval appears in the 
numerator since we are working with unitary logarithmic 
frequency interval. 
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Now we are in the position to calculate relative fluc- 
tuations SF/F of any GWB produced by astrophysical 
sources associated with galaxies within a given solid angle 
< 1. From Eqs. (|) and (|l4|) we derive 




Here the redshift z + corresponds to the beginning of star 
formation in the Universe. In our calculations we assumed 
z* = 10 (in fact, the exact value is of minor importance 
since sources at small and moderate redshifts mostly con- 
tribute to the flux). 

For stochastic GWB the flux per unit logarithmic fre- 
quency interval is related to the dimensionless strain am- 
plitude h as 




so at a given frequency fluctuations in the strain amplitude 
8h/h=(l/2)(5F/F). 

Note that the value of fluctuations depends on the spe- 
cific source type only through the dependence of energy 
carried away by gravitational waves on frequency (index 
a for the power-law dependence) and the dimensionless 

change of event rate with redshift (j^Pj ■ The depen- 
dence on frequency can appear only through the change 
in spectral index a with frequency. To ensure that the 
spectral shape has no effect on the relative fluctuations, 
observations should be performed at a frequency which 
is sufficiently (by an order of magnitude) lower than the 
high-frequency cut-off of the comoving power-law spec- 
trum of the background. 

2.3. Specific examples 

As specific examples we consider GWB produced by ex- 
tragalacting merging white dwarfs and rapidly rotating 
hot neutron stars. In the first case the energy carried 
away by gravitational radiation is the orbital energy of 
binary white dwarfs whose evolution is driven by gravi- 
tation wave emission. In the LISA sensitivity frequency 
band 1CT 4 - lfT 1 Hz we have E(f) ~ / 2 / 3 , a = 2/3. A 
white dwarf binary system reaches the orbital frequency of 
the pre-merging stage 10~ 2 Hz typically ~ 10 9 years after 
the formation, so the dependence of the rate of events on 
redshift goes smoother than star formation rate history 
SFR(z) in the Universe. The isotropic stochastic GWB 
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Fig. 1. Relative rms fluctuation of the gravitational wave 
flux from unresolved extragalactic astrophysical sources at 
different angular scales = n/q. Shown are fluctuations of 
the backgrounds produced by hot rapidly rotating neutron 
stars at frequency 100 Hz (upper curves) and by coalescing 
binary white dwarfs at frequency 10~ 2 Hz (lower curves). 
Calculations for flat FRW universes without cosmological 
constant (fl^ = 0) and with Hl — 0.7. 

produced by extragalactic binary white dwarfs with ac- 
count of the star formation rate evolution in the Universe 
was calculated by Kosenko & Postnov (1998). At the fre- 
quency 0.01 Hz its level is h w d ~ 5 x 10~ 21 . 

For hot young neutron stars the emission of gravita- 
tional waves can be driven by r-mode instability (Lind- 
blom ct al. 1998, Owen et al. 1999). The energy carried 
away by gravitational waves is the rotational energy of the 
neutron star E r and E r (f) ~ f 2 (i.e. a — 2) within the 
frequency band ~ 120 Hz < / < (2/37r)f2/f, where £Ik is 
the Kepler frequency at which mass shedding at the stellar 
equator makes the star unstable. The isotropic stochastic 
GWB produced by hot NS was studied by Ferrari et al. 
(1999b). Hot neutron stars are formed in the core collapse 
supernova explosion events in the end of evolution of mas- 
sive stars, so no deviation from star formation dependence 
on redshift for their rate is expected. The isotropic back- 
ground level at 100 Hz was estimated to be h ~ 3 x 10~ 25 . 

The relative fluctuations of these backgrounds at an- 
gular scales = n/q are shown in Fig. [I] as a function of 
wave numbers q. It is seen that at < 1 degree (q > 100), 
where the (quasi)linear regime of the galaxy clustering 
evolution with redshift is expected so that we can use 
ty(z) = (1 + z) _1 , the fluctuations amount to 5% for co- 
alescing white dwarfs and 10-20% for hot neutron stars. 
At larger scales (smaller q) the relative fluctuations slowly 
decrease. The calculations were performed for two cosmo- 
logical models: a flat FRW universe with zero cosmological 
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constant and with fl\ = 0.7. In the case of hot neutron 
stars the level of fluctuations increases in the A-dominated 
Universe at all scales considered, while in the case of coa- 
lescing white dwarfs the curves for two cosmological mod- 
els intersect at q ~ 110 (angular scales 9 ~ 1.5 degrees). 

Note that at smaller angular scales (larger q) the non- 
linear evolution of LSS should be taken into account, so 
the present calculations cannot be applied. 

3. Conclusion 

We have calculated the expected level of angular fluctua- 
tions of cosmic stochastic gravitational wave backgrounds 
produced by populations of astrophysical sources associ- 
ated with galaxies. The dependence of the source forma- 
tion rate on redshift is taken into account using global 
star formation history in the Universe. The relative rms 
fluctuations of the GW flux does not depend on the spe- 
cific source formation rate per galaxy, only on its evolution 
with redshift. The level of relative flux fluctuations at an- 
gular scales 9 < 1 degree is found to be ~ 5% for coalescing 
white dwarfs with insignificant dependence on the cosmo- 
logical model assumed, and 10—25% for hot neutron stars, 
with stronger fluctuations for flat cosmological constant 
dominated Universe. Angular dependence of these fluc- 
tuations is a distinctive feature of such GW backgrounds 
and can be used to discriminate between astrophysical and 
relic cosmological stochastic noises in future experiments. 
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